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SECTION  1 

STRESS  RATES  IN  CONTINUUM  MECHANICS  AND  COMPUTER  CODES 


The  necessity  of  paying  careful  attention  to  stress  rates  in  the 
theory  of  mechanics  of  a  continuous  medium  is  illustrated  by  a  simple 
example:  Suppose  a  body  undergoes  rigid  body  motion  (rotation  and 
translation)  without  stretching,  i.e., 


(1) 


where  x  =  position,  v  =  velocity  and  D  is  the  rate-of-deformation 
tensor.  Then,  at  every  point  in  the  body,  in  general,  the  time 
derivative  of  the  stress  tensor  is  constant  only  in  a  frame  of  refer¬ 
ence  attached  to  the  body,  i.e.,  a  rotating  frame.  In  other  words, 


(2) 


In  this  case  of  rigid  body  motion,  the  quantity 

Jj  :  J  -  1  J  ♦  J  I 

actually  is  =  0;  here  W  is  the  spin  tensor: 


(3) 
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an  infinite  family  of  stress  rat 
objectivity,  or  frame  indifference. 
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the  so-called  R i v I i n-Er icksen  tensors;  examples  of  non-frame- i nd i f- 
ferent  tensor  quantities  are  W,  v  and  a.  For  a  discussion  of  frame- 
indifference,  see  [1].  The  objective  stress  rate  typically  appears  in 
a  constitutive  relation  of  the  form 


°ij  =  Ci jk£  (°k£  "  DkJ 


where  C  is  the  elastic  constitutive  tensor  and  the  plastic  strain 
rate.  ~ 

Stress  rates  which  are  closely  related  to  the  Jaumann  rate  and 
which  have  appeared  in  the  literature  are  the  Cotter-Ri v I i n  rate: 


gCR  =  gj  ♦  g  g  ♦  g  g  , 


the  Oldroyd  rate: 


=  gj  -  (D  a  *  g  Q) 


and  the  Truesdel I  rate: 


0  0 

Oj  =  Oq  *  a  tr  D 


These  rates  are  defined  in  terms  of  tensor  quantities  evaluated  in  the 
current  material  configuration  (one  point  tensors),  and  are  easy  to 
incorporate  in  either  an  Eulerian  or  Lagrangian  computer  code.  An 
objective  stress  rate  which  is  gaining  popularity  in  Lagrangian  codes 
is  the  Green-Naghdi  rate: 


-GN 


where 


n  =  ft  rt 


and  R  is  the  proper-orthogonal  rotation  tensor  in  the  polar  decomposi¬ 
tion  of  the  deformation  gradient: 


2 


Here  X  is  the  position  of  the  particle,  currently  at  x,  in  the 
undeformed  body;  _F  is  a  so-called  two  point  tensor.  It  is  easily 

o 

seen  that  the  evaluation  of  requires  "keeping  books"  on  the 

original  position  of  each  material  point  in  the  continuum,  and  while 
feasible  in  Lagrangian  codes  would  be  considerably  more  complicated  in 
an  Eulerian  context. 

In  [2]  it  is  shown  that  of  the  Jaumann,  Oldroyd,  Cotter-Ri v I i n  and 
Truesdel I  rates,  only  the  Jaumann  rate  gives  a  linear  stress-strain 
response  for  a  hypoelastic  body  of  grade  zero  undergoing  simple  exten¬ 
sion,  i  e. ,  for 

T]  0  0 

J)  =  <M  0  (13) 

0  0  e 

where  rj  and  e  are  functions  of  time  (see  Figure  1).  In  addition,  it 
is  shown  in  [3]  that  of  the  Jaumann,  Oldroyd,  Cotter-R i v I i n  and 
Truesdel I  stress  rates,  only  the  Jaumann  rate  has  the  property  that 
when  it  vanishes  the  stress  invariants  are  stationary. 

However,  it  is  shown  in  [4]  that  the  Jaumann  rate  suffers  a 
serious  deficiency:  consider  a  linear  elastic  medium  undergoing  pure 
rectilinear  shear  (plane  strain),  with 

0  a/2  O' 

D  =  a/2  0  0  (14) 

=  0  0  0 


'  0  a/2  0 

W  =  -a/2  0  0 


and  obeying  the  constitutive  relation 


The  sinusoidal  variation  of  stress  is  clearly  physically  unreasonable. 

O 

On  the  other  hand,  it  is  shown  in  [4]  that  the  use  of  £qn  instead  of 

O 

ffj  in  (16)  gives  a  monotonic  behavior  of  stress  for  this  problem 
(Figure  2).  This  fact  has  led  to  the  implementation  of  the  Green- 
Naghdi  rate  in  Lagrangian  codes  such  as  NIKE2D,  DYNA2D  and  HONDO. 


£j  =  X  I  tr  D  ♦  2p  D  (16) 

where  X  and  /i  are  the  Lame'  constants  and  1^  is  the  identity  tensor, 
with  stress  and  strain  equal  zero  initially,  and  a  =  constant.  The 
governing  equations  are  then 


"l2  ~  0 


;i2  -  I  4 


"22  ‘  "llj 


"22  *  a  "12  *  0 


"33  *  0 


which  have  the  solution  (for  the  above  initial  conditions) 
o.0  =  p  sin  at 


a..  =  /i(l  -  cos  at) 


ff22  =  •  °ll 


0  33  =  ° 


(17a) 


(17b) 


(17c) 


(17d) 


(18a) 


(18b) 


(18c) 


(18d) 
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Figure  2.  A  comparison  between  the  shear  stress  on  a  hypo-elastic  material 
in  a  simple,  rectilinear  shear  using  the  Green-Naghdi  and  the 
Jaumann  stress  rates. 


Consider  the  above  model  problem  of  rectilinear  shear  with  a 
stress  rate  of  the  form 

j  =  !  -  #  f  ♦  a  1  ♦  a(s  D  +  D  ?)  *  P  Z  tr  D  (19) 


where  a  and  p  are  constants.  Note  that  the  Cotter-Ri v I i n,  Oldroyd  and 
Truesdell  rates  are  included  in  this  form.  For  this  particular 
problem  the  governing  equations  are 


an  ♦  (a  -  1)  a  a12  =  0 

*12  *  2 [ 4  *11+  (a  _1^a  *22]  =  ^ 

a22  ♦  (a  ♦  l)a  a12  =  0 


*33  =  0 


(20a) 

(20b) 

(20c) 

(20d) 


o\o  =  («-!)(«♦  l)a  0. 


For  the  case  (a  -  1) (a  *  1)  =0  the  solution  is 


12  3  *  a  ‘ 

(22a) 

0 

II 

H 

(a  =  1) 

(22b) 

2  .2 
=  a  fit 

(a  =  -1) 

2  .2 

22  -  -• 

(a  =  1) 

(22c) 

=  0 

(a  =  -1) 

m 


"*:**»>*«!*^ 


M 


where  is  the  plastic  strain  rate: 

D  =  D®  ♦  Dp 


where  De  is  the  elastic  strain  rate,  and 


with  yield  function 


=  X  I  tr  D®  ♦  2p  De 


f  =  1  S .  .  S..  -  \  Y2 
2  ij  ij  3 


where  Y  is  the  (constant)  stress  difference  in  compression  for  a 
von  Mises  material  and  S  is  the  stress  deviator: 


Si 


S  ■  *  “  O''  *  M  5  •  •  ffl  I 

i j  i  j  3  i  j  kk 


m 


This  is  the  von  Mises  failure  surface  and  the  Prandtl -Reuss  flow  rule. 
In  the  (initial)  elastic  solution,  J2'  =  1/2  Sjj  Sjj  reaches  a  maximum 
at  dimensionless  time  at  =  f,  and  we  suppose  plastic  failure  occurs 
before  this  point.  Subsequently,  as  long  as  the  material  is  on  the 
yield  surface,  the  governing  equations  are 


?! 

« 


“ll  -  a  "l2  =  *  S11 

512  *  2  aKl  -  "22)  *  "a  -  ^  A  S1 
*22  +  3  *12  =  "  2/i  A  S22 


(30a) 


(30b) 


(30c) 


a33  =  A  (Sn  S22^ 
Adding  (30a),  (30c)  and  (30d) ,  we  obtain 

*11  *  *22  *  *33  =  tr  * 


(30d) 


and  since  at  tne  eno  ot  tne  elastic  solution,  tr  a  ~  u,  we  nave  also 


for  the  plastic  solution  tr  a  =  0,  so  that  a  =  S. 

Now 

J2  =  Sll2  *  S222  *  S122  +  S11  S22 

(32) 

=  i  Y 2 

3 

so  that 

2SU§11  ♦  2S22§22  ♦  2S12§12  ♦  Su§22  ♦  S22Sn  =  0 

(33) 

which  gives 

-  4/i  A  (Sn2  ♦  S222  ♦  S122  ♦  Sn  S22]  ♦  2/i  a  S12  =  0 

(34) 

so  that 

A  =  3a  S12/2Y2 

(35) 

Also, 

^11  4  $22  =  -  2(i  A  (Sn  ♦  S22) 

(36) 

■  ^  (S11  *  S22>  S!2 

Since  Sji  =  -  S22  at  the  end  of  the  elastic  solution,  (36)  is  satis 
fied  thereafter  if  S11  =  -  Soo.  so  that  the  aovernina  eauations  in  th 


so  that  for  a  steady  state  to  exist,  we  must  have 


<  31'2 


Y  <  3 


The  character  of  the  solution  is  as  follows:  if:  Y  >  12*^2  fi  the 
slastic  solution  (18)  persists  forever.  If  Y  <  12*^2  fi,  the  elastic 


solution  holds  until  dimensionless  time 


-1  Y 

at  =  cos  1  -  — ~ 
o  ~  z 

6/* 


at  which  time 


S11^V  =  6/i 


S12(V  =  Y[(l  -  Y2/12^2)/3]1/2 


(41a) 


(41b) 


If  Y  <  31/2  fi,  the  solution  asymptotes  to  (38);  otherwise  the  solution 
oscillates  between  the  elastic  and  plastic  flow  regimes.  This  has 
been  verified  by  numerical  solution  of  (37).  Figure  3  shows  the  solu¬ 
tion  for  the  case  Y  =  fi  =  100  kb;  here  S^C®) /S^C^o)  =  0-85.  In  the 
case  Y  =  /i/10  (see  Figure  4),  S^2(*)/Si2(^o)  =  0-9987.  In  all  cases 
where  a  steady  state  exists,  Sjj(»)  =  2Sjj(t0). 

These  results  indicate  that  for  real  materials  with  Y  <<  fi,  the 
Jaumann  stress  rate  should  be  adequate. 


mmmt 


Figure  3.  Elastic/plastic  solution  for  Y  =  u  :  Stress  vs  Ai 
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ATTN:  DGRINE 

2  CYS  ATTN:  D  H  BROWNELL  JR 
ATTN:  LIBRARY 


CUSHING  ASSOCIATES,  INC 
ATTN:  V  CUSHING 


SCIENCE  &  ENGRG  ASSOC.  INC 
ATTN:  B  CHAMBERS 


DEVELCO.  INC 

ATTN:  LRORDEN 


SCIENCE  &  ENGRG  ASSOCIATES.  INC 
ATTN:  J  STOCKTON 


GENERAL  RESEARCH  CORP 
ATTN:  E  STEELE 
ATTN:  R  PARISSE 

GEO  CENTERS,  INC 

ATTN:  HLINNERUD 
ATTN:  L  ISAACSON 

H-TECH  LABS,  INC 

ATTN:  BHARTENBAUM 

I  IT  RESEARCH  INSTITUTE 

ATTN:  DOCUMENTS  LIBRARY 

KAMAN  SCIENCES  CORP 
ATTN:  LIBRARY 

KAMAN  TEMPO 

ATTN:  DASIAC 


SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  K  SITES 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  TECHNICAL  LIBRARY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  WLAYSON 

SOUTHWEST  RESEARCH  INSTITUTE 
ATTN:  A  WENZEL 

SRI  INTERNATIONAL 

ATTN:  DKEOUGH 
ATTN:  PDECARLI 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  D  ORMOND 
ATTN:  F  LEOPARD 


KAMAN  TEMPO 

ATTN:  DASIAC 


TERRA  TEK,  INC 

ATTN:  S  GREEN 


MERRITT  CASES,  INC  TRW  ELECTRONICS  &  DEFENSE  SECTOR 

ATTN:  J  MERRITT  2  CYS  ATTN:  N  LIPNER 

ATTN:  LIBRARY  ATTN:  TECH  INFO  CTR  DOC  ACQ 


MITRE  CORP  TRW  ELECTRONICS  &  DEFENSE  SECTOR 

ATTN:  J  FREEDMAN  ATTN:  E  WONG 

ATTN:  P  DAI 

PACIFIC-SIERRA  RESEARCH  CORP 

ATTN:  H  BRODE,  CHAIRMAN  SAGE 


DEPT  OF  DEFENSE  CONTRACTORS  (CONTINUED) 

WASHINGTON  STATE  UNIVERSITY 
2  CYS  ATTN:  PHYSICS  DEPT  Y  M  GUPTA 

WEIDLINGER  ASSOC,  CONSULTING  ENGRG 
ATTN:  TDEEVY 


WEIDLINGER  ASSOC,  CONSULTING  ENGRG 
ATTN:  M  BARON 

WEIDLINGER  ASSOC.  CONSULTING  ENGRG 
ATTN:  J  ISENBERG 
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